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Abstract 

We investigated the variability of the climate-growth relationship of Aleppo pine across its distribution range in the 
Mediterranean Basin. We constructed a network of tree-ring index chronologies from 63 sites across the region. Correlation 
function analysis identified the relationships of tree-ring index to climate factors for each site. We also estimated the 
dominant climatic gradients of the region using principal component analysis of monthly, seasonal, and annual mean 
temperature and total precipitation from 1,068 climatic gridpoints. Variation in ring width index was primarily related to 
precipitation and secondarily to temperature. However, we found that the dendroclimatic relationship depended on the 
position of the site along the climatic gradient. In the southern part of the distribution range, where temperature was 
generally higher and precipitation lower than the regional average, reduced growth was also associated with warm and dry 
conditions. In the northern part, where the average temperature was lower and the precipitation more abundant than the 
regional average, reduced growth was associated with cool conditions. Thus, our study highlights the substantial plasticity 
of Aleppo pine in response to different climatic conditions. These results do not resolve the source of response variability as 
being due to either genetic variation in provenance, to phenotypic plasticity, or a combination of factors. However, as 
current growth responses to inter-annual climate variability vary spatially across existing climate gradients, future climate- 
growth relationships will also likely be determined by differential adaptation and/or acclimation responses to spatial 
climatic variation. The contribution of local adaptation and/or phenotypic plasticity across populations to the persistence of 
species under global warming could be decisive for prediction of climate change impacts across populations. In this sense, a 
more complex forest dynamics modeling approach that includes the contribution of genetic variation and phenotypic 
plasticity can improve the reliability of the ecological inferences derived from the climate-growth relationships. 
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Introduction 

Climate strongly influences the geographical distribution of 
plant species in terrestrial ecosystems [1,2]. As local temperature 
and water availability change with broader changes in climate, 
biodiversity, species distribution, and growth rates are expected to 
change accordingly [3,4,5]. However, the impact of climate 
change may vary gready depending on the ability of species to 
acclimate or to adapt to future climate conditions [6]. 

Plant species can adjust to new environmental conditions 
through local adaptation or phenotypic plasticity [6]. Local 
adaptation may imply genetic differentiation among populations 
as a consequence of differential selection pressures and/or 
population isolation [7]. Phenotypic plasticity, defined as the 
range of phenotypes that a single genotype can express as a 
function of its environment, can also be a crucial factor for plant 



response to rapid climate change [6,8]. Heterogeneous responses 
of species to climate variability across their range are directly 
connected to phenomena of local adaptation and phenotypic 
plasticity and are the basis of potential adaptability to future 
climate conditions. Detailed knowledge of the relationship 
between climate and growth across the range of distribution of 
species is essential to predict and mitigate the effects of climate 
change. 

Climate change is especially rapid and extreme in the 
Mediterranean basin [9,10,11]. However, due to the transitional 
nature of the Mediterranean climate (which ranges from near- 
desert to temperate regimes), Mediterranean areas contain a great 
variety of natural conditions which will likely affect the impact of 
climate change [12,13,14]. 

In the widely diverse set of climatic conditions contained in the 
Mediterranean basin, pines are by far the most widespread genus 
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covering approximately 5% of the total land area and 25% of the 
forested area. The most common species is the Aleppo pine (Pinus 
halepensis Mill.) occupying large areas of the western Mediterra- 
nean as well as occurring in the eastern Mediterranean [15,16]. 
The ecology and biogeography of Aleppo pine is well-researched 
[17,18]. The dendroclimatology of Aleppo pine has been 
investigated at individual locations within the Mediterranean 
basin [19,20,21,22,23,24,25,26,27,28,29] confirming that Aleppo 
pine is sensitive to climate variation. However, a comprehensive 
dendroclimatic analysis across the distribution range of Aleppo 
pine is lacking. 

In this paper, we investigate the variability of the climate-growth 
relationships of Aleppo pine using a dendroclimatic network 
composed of high-resolution climatic and dendrochronology 
databases comprising the complete climatic gradient across the 
Mediterranean Basin. Our objectives were to a) characterize the 
climatic variability within the range of Aleppo pine; b) relate high 
frequency variation in tree-ring width to monthly, seasonal, and 
annual temperature and precipitation; c) determine the depen- 
dence of the dendroclimatic relationships to underlying climatic 
patterns. 

Methods 

Ethics Statement 

All sampling sites were located in public forests. The necessary 
permits for field sampling were issued by the forestry authorities at 
regional and local level. The locations were not protected areas, 
and the field studies did not involve endangered or protected 
species. 

Locations and climatic data 

The distribution map for Aleppo pine (Fig. 1) based on an 
earlier map [30] was obtained from the European Forest Genetic 
Resources Programme (EUFORGEN http://www.euforgen.org/ 
distribution_maps.html). Climatic time series data for 1,068 spatial 
grid points within the distribution range for the period 1901-2000 
were obtained from the Climatic Research Unit (CRU) of the 
University of East Anglia. For the western part (Spain, Algeria, 
France, Italy, Slovenia, Greece, and Turkey), monthly, seasonal, 
and annual mean temperatures and total precipitation were 
provided by the CRU TS 1.2 dataset (10-minute resolution) [31]. 
For the portion of the distribution located east of the coverage of 
CRU TS 1.2, we used the CRU TS 2.1 dataset (0.5° grid 
resolution) [32]. To characterize the spatial variability of the 
climate across the distribution area of Aleppo pine, while 
minimizing the effect of outlier values, we determined the central 
95% range for each variable, defined as the values occurring 
between the 2.5 and 97.5 ranked percentiles. 

The dendrochronological network consisted of newly collected 
and archived tree-ring width series from 63 Aleppo pine sites 
(Table 1), between 32.23° and 45.67°N latitude, 1.41°W and 
36.17°E longitude, and 15 to 1650 m a.s.l. (Fig. 1). To determine 
the adequacy of the network to describe the range of dendrocli- 
matic responses, climatic data for the CRU gridpoint closest to 
each site were compared with those of the Aleppo pine distribution 
as a whole. 

A T-mode Principal Component Analysis (PCA) [33] was used 
to summarize the spatial variability of the mean climate conditions 
across the study area. PCA is a data reduction technique that 
transforms a large group of variables into a new smaller set of 
variables called principal components (PC), which are linear 
combinations of the original variables. These PC were calculated 
on the correlation matrix of mean monthly temperatures and 



mean total monthly precipitation for the period 1901-2000. The 
components were rotated (Varimax) to redistribute the explained 
variance and to obtain more stable and robust spatial patterns 
[34]. Selected PC, guided by Kaiser's Rule (eigenvalues >1; [35]), 
described climatic gradients across the Mediterranean basin were 
used in subsequent analysis. 

Tree-ring chronology construction 

The tree-ring width series were derived from a total of 1634 
increment cores from 818 trees (Table 1). Tree selection, core 
collection, processing, and ring-width measurement were con- 
ducted using standard techniques [36,37]. Ring-width index 
chronologies for each site in the network were constructed using 
ARSTAN software version 6.05P [38] to retain the high-frequency 
variation that would be most sensitive to the climate variables to be 
tested. First, a negative exponential or linear regression function 
was fitted to each ring-width series. The residuals from this first 
detrending were then fitted with a cubic smoothing spline function 
(50% frequency cut-off of 30 years). The residuals from the second 
detrending were then fitted with an autoregressive model to reduce 
the autocorrelation in the model residuals which were then 
averaged to construct the ring-width residual index chronology 
(chronology RES in ARSTAN) for each site using the biweight 
robust mean function. An expressed population signal (EPS) of 
0.85 was used as a threshold to determine the reliable part of each 
index chronology to use in subsequent analysis [39] . 

Dendroclimatic relationships across climate gradients 

For our first-stage tests, the statistical relationship between 
monthly, seasonal and annual climate series and annual tree-ring 
chronologies was assessed individually for each of the 63 sites by 
correlation function (CF) analysis using the program Dendro- 
Clim2002 [40] . The RES chronology was the dependent variable 
while the regressors were the monthly, seasonal and annual mean 
temperatures and the total precipitation for each 16-month 
biological year (from the previous September to the current 
December). 

The second-stage test consisted of multiple regression analysis 
(forward selection) using the CF coefficients calculated in the first 
stage as dependent variables and the significant principal 
components (PC) as independent variables that characterize mean 
climate conditions for each of 63 study sites. In this test, 
dependence of CF on PC indicates that the correlation of growth 
to annual, seasonal, and monthly climate variability depended on 
the underlying climatic conditions. Then, if dendroclimatic 
relations were found to be climate-dependent, the variation in 
the climate-growth relationship may be quantified and predicted 
for the whole distribution area of species by resolving obtained 
regression models on all 1,068 points for which PC values were 
calculated. 

Results 

Climate characteristics within the distribution area of 
Aleppo pine 

The climatic conditions varied widely from near desert to moist 
temperate. The central 95% of the range (bounded by the 0.025 
and 0.975 percentiles) of mean annual temperatures ranged from 
10.2°-17.9°C. Winter temperature (average of 7.0°C) ranged from 
2.9-1 1.5°C and summer temperature (average of 22.3°C) ranged 
from 18.2-25.9°C (Fig. 2). The central 95% range of annual total 
precipitation (average of 570 mm) across the geographic distribu- 
tion ranged from 328-957 mm. The seasonal distribution of 
precipitation varied greatly with winter precipitation comprising 



PLOS ONE | www.plosone.org 



2 



December 2013 | Volume 8 | Issue 12 | e83550 



Dendroclimatic Variability in Aleppo Pine 




Figure 1. Map of the Mediterranean Basin and the distribution of Aleppo pine (A); details in B, C, D (bottom). 

doi:10.1371/journal.pone.0083550.g001 



about 3 1 % of the total annual precipitation. Among sites, winter 
precipitation contributed 15-66% of the annual total precipita- 
tion. Spring precipitation contributed an average of about 27% of 
the total annual precipitation with sites ranging from about 12— 
36%. Autumn precipitation contributed an average of about 30% 
with a range of about 13-43%. Summer was generally the dry 
season with an average annual contribution of about 1 1 % of the 
total precipitation but ranged across sites from less than 0.1% to 
about 29% of the total annual precipitation (Fig. 2). 

The climate variation across the distribution area of Aleppo 
pine can be summarized by four significant principal components 
(PC) which together explain about 95% of the variation. PCI 
(containing 59% of the variation) represented the main climatic 
gradient, from cold and wet conditions to warm and dry (Fig. 3). 
PC2 (about 24%) related to variation in precipitation (Fig. 3). PC3 
(about 8%) represented variation in the seasonal distribution of 
precipitation from sites where winter is the main precipitation 
season to sites where the contribution of summer and autumn 
precipitation to the total annual precipitation increases (Fig. 3). 
PC4 (about 3%) represented variation in temperature between 
winter and summer (Fig. 3). 

Dendroclimatic relationships in the Aleppo pine tree-ring 
network 

The basic statistics for 63 local chronologies are shown in 
Table 1. Mean correlation between trees and mean sensitivity of 
chronologies varied from 0.30 to 0.83 (average of 0.59) and from 



0.18 to 0.65 (average of 0.34) respectively, indicating the potential 
for a common climate signal of varying strength among sites. 

The CF analysis showed that ring-width index was primarily 
related to precipitation and secondarily to temperature. However, 
a high variability in the relationship between climate and growth 
was observed across the network (Fig. 4). 

Tree-ring index was significantly related to annual total 
precipitation for all but two chronologies (FAR in France and 
AMA in Greece). The correlation of tree-ring index to winter 
temperature was usually significant and positive in colder 
locations, indicating reduced growth in years with cold winters. 
In contrast, at three network sites, where conditions are warmer, 
significant negative correlations indicated reduced growth in years 
with warm winters. 

Variation in dendroclimatic relationships of Aleppo pine 
across its distribution 

The significant regression equations of CF with the four PCs 
indicated that dendroclimatic relationships depended on the site 
position along the climatic gradient of the network (Table 2). This 
dependence produced a clear and coherent spatial pattern of tree 
growth in the Mediterranean Basin (Fig. 5). The CFs of annual 
temperature were inversely related to PC 1 indicating that reduced 
growth associated with warm conditions occurred mainly in the 
southern part of the species distribution area where temperature 
was generally higher and precipitation lower than the regional 
average. Growth limitations due to cold conditions occurred in the 
northern part of the network where the average temperature was 
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Table 1. Study sites. 





N° 


Site 
code 


Site 


Country 


Lat 


Long 


Elevation 


Trees 


Samples 


Tree- 
rings 


EPS>0.85 


MS 
(res) 


MC 
(res) 


1 


BOU 


Bout10(***) 


Algeria 


35.65 


5.5 


1200 


11 


26 


2463 


1893-1988 


0.3 


0.43 


2 


CHE 


Chelia 9(***) 


Algeria 


35.25 


6.85 


1650 


11 


32 


4330 


1819-1988 


0.25 


0.68 


3 


BEZ 


Bezez 8(**») 


Algeria 


35.2 


6.9 


1300 


12 


34 


3558 


1864-1988 


0.4 


0.67 


4 


TZI 


Tizi Tmellet 6(**») 


Algeria 


35.13 


6.8 


1450 


10 


24 


3264 


1831-1988 


0.47 


0.74 


5 


RAS 


Ras fourar 7(***) 


Algeria 


35.05 


6.7 


1200 


9 


27 


3342 


1835-1987 


0.47 


0.72 


6 


SAH 


Sahari 3(***) 


Algeria 


35.05 


3.5 


1060 


11 


32 


3843 


1849-1989 


0.3 


0.6 


7 


MES 


Messaad 5(***) 


Algeria 


35 


4.1 


1100 


11 


31 


3354 


1797-1989 


0.41 


0.75 


8 


SAF 


Safai 4(***) 


Algeria 


34.9 


3.9 


1100 


11 


28 


2925 


1851-1989 


0.4 


0.72 


9 


SEN 


SenalbaNf**) 


Algeria 


34.72 


3.12 


1410 


11 


33 


3322 


1858-1988 


0.31 


0.68 


10 


SEW 


SenalbaW(***) 


Algeria 


34.65 


3.08 


1350 


11 


33 


3451 


1859-1990 


0.37 


0.74 


11 


NYO 


NyonsH 


France 


44.36 


5.16 


350 


8 


14 


1171 


1906-1993 


0.23 


0.6 


12 


MEE 


Les Mees(*) 


France 


44.03 


5.99 


600 


9 


12 


941 


1912-1994 


0.29 


0.58 


13 


PER 


Peruis(*) 


France 


44.03 


5.91 


600 


8 


13 


1001 


1924-1992 


0.18 


0.43 


14 


ROB 


RobionP) 


France 


43.84 


5.12 


300 


11 


20 


1663 


1905-1994 


0.31 


0.57 


15 


CVI 


Vitrolles(*) 


France 


43.82 


5.58 


600 


9 


12 


919 


1918-1993 


0.28 


0.51 


16 


MOU 


Moustier Sainte-MarieH 


France 


43.82 


6.21 


750 


10 


12 


1062 


1900-1994 


0.32 


0.57 


17 


FOO 


MerindolH 


France 


43.78 


5.18 


330 


9 


13 


1398 


1881-1994 


0.35 


0.69 


18 


RDA 


Roque d'antheron(*) 


France 


43.71 


5.33 


170 


13 


19 


1645 


1896-1994 


0.27 


0.59 


19 


RIA 


RiansH 


France 


43.6 


5.75 


450 


10 


14 


1179 


1903-1994 


0.32 


0.63 


20 


ROU 


Rousset(*) 


France 


43.5 


5.62 


200 


6 


6 


426 


1920-1990 


0.31 


0.71 


21 


ROG 


Rognac(*) 


France 


43.48 


5.26 


190 


7 


15 


1346 


1903-1994 


0.34 


0.56 


22 


PNM 


Les Pennes Mirabeau(*) 


France 


43.4 


5.29 


150 


7 


14 


1050 


1910-1993 


0.34 


0.6 


23 


AUR 


Auriol(*) 


France 


43.36 


5.67 


300 


8 


11 


1040 


1901-1993 


0.28 


0.51 


24 


MAR 


Marseille(**) 


France 


43.22 


5.27 


200 


7 


21 


2664 


1833-1973 


0.42 


0.64 


25 


CIO 


La Ciotatf) 


France 


43.21 


5.58 


200 


12 


17 


1517 


1899-1994 


0.31 


0.5 


26 


FAR 


Toulon(*) 


France 


43.15 


5.94 


500 


9 


15 


1353 


1904-1993 


0.23 


0.45 


27 


KAS 


Kassandra 


Greece 


39.98 


23.47 


136 


10 


30 


1369 


1937-1989 


0.27 


0.71 


28 


KSY 


Athensl 


Greece 


38.05 


23.82 


260 


10 


17 


1777 


1910-2003 


0.31 


0.48 


29 


DBG 


Athens2 


Greece 


38 


23.62 


180 


23 


37 


3183 


1909-2001 


0.35 


0.37 


30 


AMA 


Amaliada 


Greece 


37.9 


21.4 


88 


12 


36 


1400 


1948-1989 


0.21 


0.64 


31 


CAR 


Carmel Mountain(*****) 


Israel 


32.67 


35 


450 


8 


9 


613 


1959-1997 


0.27 


0.43 


32 


TER 


Terni 


Italy 


42.62 


12.65 


470 


33 


33 


1985 


1940-2003 


0.21 


0.3 


33 


DIB 


Dibeen(*****) 


Jordan 


32.23 


35.82 


867 


14 


14 


819 


1936-1993 


0.29 


0.4 


34 


SLO 


Krkavce-Dekani 


Slovenia 


45.67 


13.7 


50 


25 


49 


2310 


1949-2004 


0.23 


0.38 


35 


AYE 


Ayerbe (Biel) 


Spain 


42.32 


-0.84 


924 


19 


33 


1377 


1962-2006 


0.23 


0.45 


36 


EBA 


Ejea-Bardenas 


Spain 


42.17 


-1.4 


365 


7 


7 


499 


1941-2003 


0.32 


0.56 


37 


GRA 


El Grado 


Spain 


42.16 


0.2 


168 


15 


30 


1506 


1950-2006 


0.24 


0.56 


38 


EST 


Estopinan del Castillo 


Spain 


41.97 


0.61 


502 


14 


27 


1087 


1965-2006 


0.38 


0.62 


39 


VLL 


Villanueva de Gallego 


Spain 


41.88 


-0.91 


452 


15 


29 


2692 


1898-2006 


0.4 


0.58 


40 


CAP 


Alcubierre (San Caprasio) 


Spain 


41.75 


-0.5 


738 


14 


27 


1516 


1942-2006 


0.29 


0.63 


41 


FRA 


Fraga 


Spain 


41.47 


0.32 


340 


14 


29 


3858 


1857-2006 


0.43 


0.62 


42 


CAS 


Caspe 


Spain 


41.29 


0.07 


166 


15 


28 


3539 


1856-2007 


0.59 


0.68 


43 


CHI 


Chiprana 


Spain 


41.24 


-0.09 


160 


7 


9 


756 


1933-2003 


0.32 


0.45 


44 


DAR 


Daroca 


Spain 


41.14 


-1.41 


937 


14 


28 


1699 


1936-2006 


0.42 


0.83 


45 


OLI 


Oliete 


Spain 


40.99 


-0.69 


530 


15 


28 


1099 


1961-2006 


0.3 


0.5 


46 


ALL 


Alloza 


Spain 


40.98 


-0.57 


595 


17 


31 


2878 


1905-2006 


0.37 


0.69 


47 


ZOR 


Zorita 


Spain 


40.74 


-0.11 


857 


15 


30 


3004 


1884-2001 


0.4 


0.61 


48 


ORO 


Oropesa 


Spain 


40.06 


0.12 


1 


15 


30 


2022 


1929-2003 


0.32 


0.54 


49 


MON 


Montanejos 


Spain 


40.06 


-0.54 


569 


16 


31 


1223 


1956-2001 


0.39 


0.67 
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Table 1. Cont. 







N° 


Site 
code 


Site 


Country 


Lat 


Long 


Elevation 


Trees 


Samples 


Tree- 
rings 


EPS>0.85 


MS 
(res) 


MC 
(res) 


50 


GIL 


Gilpt Rant Fsnirttl 

UNCI ^Jcllll Lj|JIIIL/ 


^— — 


39 67 


— 0 35 


1 75 


14 


27 


2323 


1Q1 0-7006 

1 _7 1 U Z.UUU 


0 65 


0 7 


51 


REQ 


R pci i ip n 3 
ncuut: la 


pain 


39 47 


— i 2 


721 


1 5 


31 


3858 


1816-7003 


0 35 


0 48 


52 


JAL 


\a Id nCG 


pain 


39 1 9 


— 1 1 5 


571 


1 5 


35 


3034 


1881 -7003 

1 OO 1 ZUUJ 


0 41 


0 65 


53 


JAV 


Javea 


Spam 


38 73 


0 19 


96 


1 5 


60 


3434 


1Q97 ?nnn 

1 7Z/ — ZUUU 


0 35 


0 73 


54 


ALC 


Alcoy (Penacjuila) 




38 68 


— 0 36 


674 


1 5 


30 


3792 


1 OOj ZUU 1 


0 32 


0 51 


55 


FNT 


ttlLUy v rUML nUJdJ 


Spain 


38 67 


q 54 


1 022 


1 4 


24 


1 960 


irri 7nnft 

1 OO I — ZUUU 


0 35 


0 61 


56 


BIA 


Biar 


Spain 


38.62 


-0.77 


806 


15 


29 


1798 


1930-2001 


0.29 


0.73 


57 


MAI 


Maigmo 


Spain 


38.52 


-0.64 


845 


15 


31 


2697 


1896-2000 


0.4 


0.65 


58 


CRE 


Crevillente 


Spain 


38.29 


-0.77 


285 


12 


43 


3715 


1911-2000 


0.45 


0.49 


59 


GUA 


Guardamar 


Spain 


38.1 


-0.65 


15 


30 


77 


6047 


1916-2006 


0.36 


0.59 


60 


FUE 


Fuensanta 


Spain 


37.94 


-1.12 


138 


14 


26 


2094 


1908-2007 


0.44 


0.64 


61 


SEN 


Sierra Espuha 


Spain 


37.86 


-1.52 


846 


16 


29 


2670 


1907-2007 


0.25 


0.49 


62 


CAT 


Cartagena 


Spain 


37.61 


-1.01 


116 


15 


27 


2039 


1921-2007 


0.42 


0.61 


63 


BAB 


Bayat Bademleri(****) 


Turkey 


37.02 


30.28 


700 


15 


15 


2567 


1752-2001 


0.21 


0.61 



•DendroDB (Nicault); 
"ITRDB (Serre-Bachet); 
***DendroDB (Safar); 
****ITRDB (Schweingruber), 
*****ITRDB (Touchan), 

General information on the 63 study sites in the Aleppo pine network. LAT: Latitude ( N); LON: Longitude (°E); ELE: Elevation (m a.s.l); MS: mean sensitivity; MC: mean 
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lower and precipitation less limiting than the regional average 
(Fig. 5). During winter, the area where cold conditions were 
associated with reduced tree growth was extended to the south. In 
contrast, in spring and summer reduced tree growth associated 
with warm conditions occurred in extended areas towards the 
southern distribution limits of Aleppo pine (Fig. 5). 

Similarly, the inverse relationship of CF to PC2 suggests that the 
positive effects of annual, winter, and spring precipitation on tree 
growth significandy increased from wetter to drier sites. The 
influence of winter precipitation also increased along the wet to 
dry gradient being especially important at sites with high annual 
temperatures (Fig. 5). PC3 and PC4 were also significantly related 
to CFs suggesting that the variation in the seasonal distribution of 
precipitation and temperature across the network influenced the 
observed dendroclimatic relationships. PC 3 was positively related 
with summer and negatively with winter precipitation while PC4 
was positively related with sites with especially high summer 
temperatures. 

Discussion 

Sensitivity of Aleppo pine to climate variations across its 
distribution area 

The dendrochronological network of 63 sites comprises an 
essentially complete climatic gradient for Aleppo pine and 
provides information on how the growth of this species responds 
to climate variability across its distribution range. In the widely 
diverse set of climatic conditions, our study demonstrated that 
Aleppo pine is an excellent species for climate-growth studies. The 
mean correlation between trees in essentially all chronologies 
indicated that despite micro-environmental influences, a common 
pattern of year-to-year variation in ring width exists across all sites. 
The mean tree-ring sensitivity (the average relative difference in 



adjacent rings and a measure of climatic responsiveness) was 
generally higher than observed for other species in the Mediter- 
ranean area [41,42,43,44,45,46] highlighting the importance of 
Aleppo pine as a climate proxy for large regional areas. This could 
be especially important in the Mediterranean region where 
available climate observations are generally limited to the second 
half of the 20 th century [47], Additional sampling sites at locations 
predicted to contain strong climatic signals based on these results 
would be especially useful to improve reconstruction of past 
climate. Distribution maps of climate-growth relationships across 
the distribution range of Aleppo pine may be also especially 
valuable to predict the response of trees to changing climate 
patterns. 

High sensitivity of a species to climate variations across a wide 
geographic distribution area or in a wide range of climate 
conditions as observed here and in other dendrochronological 
studies suggests that climate change impacts may not be restricted 
to the ecotone or edge but may occur across the whole distribution 
range [48,49]. According to this, our results for Aleppo pine 
suggest that in the colder sites where low winter temperatures are 
presently limiting, future increased temperature may have a 
positive effect on tree growth. However, in the drier and warmer 
southern Mediterranean, where current positive anomalies in 
temperature have negative effects on tree growth, conditions for 
future tree growth will become less favorable as temperatures 
increase. 

The role of phenotypic plasticity in a species response to 
climate variability under global warming 

The ability of a species to respond through expression of existing 
genetic variation or phenotypic plasticity to changing environ- 
mental conditions may play a decisive role in species persistence or 
expansion under future global warming [8,50]. Phenotypic 
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Figure 2. Climate characteristics across the distribution area of Aleppo pine. Dark grey and light grey areas represent climate envelopes 
defined by the 2.5th-97.5th and the 0.5th-99.5th percentiles, respectively. Box plots show the range of mean monthly air temperatures and total 
monthly precipitation for the P. halepensis dendrochronological network of 63 locations. The central vertical lines indicate median values, boxes 
enclose the central two quartiles, whiskers indicate the 10th and the 90th percentiles and the dots represent the full range of mean climate values. 
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plasticity may be especially important as an adaptive strategy in 
trees and woody plants since the long generation time of many 
perennial species implies that the same set of genotypes needs to 
cope with year-to-year changing environmental conditions [51]. 

Generally, tree species are considered to express moderate to 
high plasticity in their responses to environmental stress [52,53]. 
Highly variable species that are able to survive in a broader range 
of environmental conditions are expected to better adjust to future 
climate conditions [7]. 

Pinus halepensis shows a wide ecological breadth and is adapted to 
a large range of environmental conditions, abiotic stressors, and 
perturbations [54]. Several recent studies highlight substantial 
phenotypic plasticity of Pinus halepensis in relation to different 
anatomical, reproductive and vegetative traits 
[55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70]. 

A substantial plasticity in the annual rhythms of cambial activity 
of Aleppo pine in response to different climatic conditions has 
been also well established. Thus, some studies suggest that the 
cambium is able to maintain activity throughout the whole year 
when and where climate conditions are favourable [71]. Under 
other conditions, cambial activity can stop for one to three months 
during winter, depending on the prevalence of low temperatures 
[72] . In addition, cambial activity in Pinus halepensis can slow down, 
or even stop, during summer drought [59,71,72,73,74] and 



resume later when moisture availability increases in autumn 
[59,71,73,74,75,76]. In other circumstances, the trees can be 
subjected to "double stress" characterized by two stops in cambial 
activity - one during winter, caused by low temperatures, and one 
during the summer triggered by high temperatures and lack of 
precipitation [73,74,77,78]. 

However, despite such evidence, plasticity in species responses 
to climate variability has been poorly explored from a dendro- 
chronological perspective and variability in dendroclimatic 
response is more often considered as due to environmental noise 
in the tree-ring signal rather than a consequence of the plastic 
character of the species. 

Geographical variation in climate conditions across the 
distribution area of a species creates differential selection pressures, 
and as a consequence population responses to climate variability 
are likely to vary among genetic provenances as well as through 
phenotypic plasticity. In this sense, our results highlight significant 
geographical variations in growth response of Aleppo pine to 
interannual climate variability suggesting substantial plasticity of 
growth in response to various climatic conditions. These results do 
not resolve the source of variability as genetic variation in 
provenance, phenotypic plasticity or some combination of these or 
other factors. However, the use of isoenzymes, terpene composi- 
tion, and other markers has shown that Aleppo pine contains little 
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Figure 3. Spatial distribution of the first four significant principal components (PC) reflecting the spatial variability of the mean 
climate conditions (left). Component loading of each PC against mean monthly, seasonal and annual temperatures and precipitation for the 
period 1901-2000 (right). 
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genetic variability [79,80] and that intrapopulation genetic 
variability is generally much higher than interpopulation variabil- 
ity [81]. These findings suggest that phenotypic plasticity has a 
dominant role in the observed variation in dendroclimatic 
response across distribution range of Aleppo pine. Then, as 
current growth responses to inter-annual climate variability vary 
spatially across existing climate gradients, we expect that future 
climate-growth relationships will also vary with position along 
future climate gradients. 

Implications 

Historical records of biological responses to climate in situ are 
perhaps the most dependable data we have to reconstruct past 
climate conditions, to predict future climate change impacts, and 
to develop realistic mitigation strategies. The tree-ring record is 
the most important and widely used source of long-term proxy 
data as they are available over a wide range of temporal and 
spatial scales. 

However, the use of tree-rings to reconstruct past climate 
patterns and to forecast of future forest growth within the context 
of global warming is based on the uniformitarian principle (UP) 



that assumes that the climate-limiting factors that controlled tree- 
ring characteristics in the past continue the present and will extend 
into the future [82]. Under UP, the modeled dendroclimatic 
relationship is viewed as stationary and consistent through time 
and the range of climatic variability. 

To the contrary, our results show that dendroclimatic relation- 
ships significantly vary across species distribution in accordance to 
underlying climate conditions. Our results also suggest that such 
variations are likely to be related to the acclimation of trees to the 
new environmental conditions through a plastic response over 
decadal time periods. This complicates the simple projection of 
current relationships into the past or the future. For example, the 
response to winter temperature change from positive to negative 
moving from northern (mesic) to southern (drier) sites, with no 
response at intermediate sites along the thermal gradient. In this 
sense, the recent warming has a positive effect at the northern limit 
and a negative effect on the southern portion of the distribution 
range. At the same time, as summer drought severity increases the 
response to precipitation disappears, as the tree enters a quiescent 
and nonreactive state. 
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Figure 4. Correlation between the site chronologies of Aleppo pine and climate; monthly, seasonal, and annual precipitation and 
temperature from September of previous year (Sep-1) to December (Dec) of the current year (top). Only significant values are shown 
(p<0.05). Box plot showing the variability of correlation coefficients across the chronology network. The central horizontal lines indicate the median 
values, boxes enclose the central two quartiles, whiskers indicate the 10th and the 90th percentiles and the dots represent the full range of 
correlation coefficients. Dark and light grey areas indicate significance levels at 99 and 95%, respectively. 
doi:10.1371/journal.pone.0083550.g004 
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Figure 5. Correlation coefficients across the distribution area of Aleppo pine calculated by using the multiple regression models 
detailed in Table 2, applied to the values of the significant principal components (PC). 
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Table 2. Variation of dendroclimatic relationships across climate gradients. 
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Accordingly, the application of UP to dendrochronology may 
be inappropriate when the dendroclimatic responses are likely to 
be inherently unstable and climate-dependent as we have shown 
for Aleppo pine in the Mediterranean basin. 

Interestingly, a number of tree-ring studies have addressed 
changes in tree sensitivity and/ or changes in the response of tree 
growth to climate in coincidence with unprecedented climate 
warming over the recent decades [83,84,85,86,87,88,89,90]. The 
causes of this lack of stability in the dendroclimatic relationship are 
still not well understood as all potential forcing factors (e.g., 
changes in climate, atmospheric C02 concentration, or nitrogen 
deposition) capable of driving this change covary and obscure the 
individual impacts [86] but a plausible explanation could be the 
acclimation of trees to the new environmental conditions through 
a plastic response over decadal time periods. Our results are 
related to a limited geographical area and a specific dataset and 
may have limited application to other species or regions. However, 
such patterns may occur elsewhere and if so, it may have serious 
implications by affecting reliability of tree-ring based climate 
reconstructions. 

In addition, climate model validation [91], species distribution 
models [7,9] and forecasts of future forest growth within the 
context of global warming [1,22,92,93,94] whether or not based 
on dendrochronology also assumes a uniform response of species 
to climate variability and change without the competitive 
advantage conferred by genetic and/or phenotypic plasticity 
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are telling us in a changing world. 
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